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Abstract 

We show that the expression for the supersymmetric partition function 
of the chiral Unitary (Laguerre) Ensemble conjectured recently by Splittorff 
and Verbaarschot [13] follows from the general expression derived recently by 
Fyodorov and Strahov [16]. 

A class of random matrices that has attracted a considerable attention recently 
[1, 2, 3, 4, 5, 6, 8, 7, 9, 10] is the so-called chiral (Gaussian) Unitary Ensemble 
(chGUE), also known as the Laguerre ensemble. The corresponding matrices are 
/ W \ 

of the form D = I - ^ I , where W stands for a complex matrix, with 1U T 

being its Hermitian conjugate. The off-diagonal block structure is characteristic for 
systems with chiral symmetry. The chiral ensemble was introduced to provide a 
background for calculating the universal part of the microscopic level density for 
the Euclidian QCD Dirac operator, see [11] and references therein. Independently 
and simultaneously it was realised that the same chiral ensemble is describing a new 
group structure associated with scattering in disordered mesoscopic wires [2]. One 
of the main objects of interest in QCD is the so-called Euclidean partition function 
used to describe a system of quarks characterized by nf flavors and quark masses 
rrif interacting with the Yang-Mills gauge fields. At the level of Random Matrix 
Theory the true partition function is replaced by the matrix integral: 

rif 

Z nf (M f ) = [VWU det{W + mfh^je-^i^) (1) 

J k=l 

where M/ = diag (jn^\ mlp^ and V(z) is a suitable potential. Here the integra- 
tion over complex W replaces the functional integral over gauge field configurations 



1 



[11]. Then the calculation of the partition function amounts to performing the 
ensemble average of the product of characteristic polynomials of iD over the prob- 
ability density P(W) oc e ~ Ar Try(w t w / )^ j n ^ e general case of non-zero topological 
charge v > the matrices W have to be chosen rectangular of size N x (N + v) 
[11]. For simplicity one may choose the probability distribution to be Gaussian as 
defined by the formula: dV(W) oc dWdW^ exp - [nTtW^W . 

The characteristic feature of the chiral ensemble is the presence of a particular 
point A = in the spectrum, also called the "hard edge" [3]. The eigenvalues of 
chiral matrices appear in pairs ±A& , k = 1,...,N. Far from the hard edge the 
statistics of eigenvalues is practically the same as for usual GUE matrices without 
chiral structure, but in the vicinity of the edge eigenvalues behave very differently. 

Let Z^[m] be the following spectral determinant (characteristic polynomial of 
it)): 

Z N [m] = det (m 2 l N + W ] w) . (2) 



and let us consider a more general (supersymmetric) 
chGUE defined as 



/C(M/, M b ) 
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This correlation 



function contains much more information on spectra of chiral matrices than the 
partition function (1) since it involves both product and ratios of the characteristic 
polynomials. Many efforts were spent on developing methods allowing one to cal- 
culate particular cases of such a general supersymmetric partition (or correlation) 
function [9, 10]. In particular, the case v — was completely solved in [15] by a 
variant of the supermatrix method [12, 10] augmented with a generalization of an 
Itzykson-Zuber type integrals [4, 5] to integration over non-compact group mani- 
folds. In the recent paper [13] Splittorff and Verbaarschot conjectured the result for 
arbitrary integer v > in the microscopic (sometimes also called " chiral" ) large- N 
limit: N — > oo such that Xbj = 2NM(,j is finite. The authors of [13] used the 
advanced version of the replica method suggested recently by Kanzieper [14]. The 
final result is given in terms of a determinant containing modified Bessel functions 
Ii(z) ("compact integrals") and their noncompact partners - Macdonald functions 
Ki(z). The goal of the present Letter is to show that the case v ^ considered 
in [13] in fact follows from a very general expression derived in the recent paper 
[16]. The demonstration of this fact also provides a natural explanation why both 
compact and non-compact integrals must appear on equal basis. 
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The eigenvalues xi, xn of the N x N positive definite matrix H = W^W are 
known to be distributed according to the Laguerre ensemble density function 



N 



V(x 1 ,...,x N ) ^ A 2 (X)Y[w u (x t ) 



(4) 



i=l 



where w v (x) = x u e~ Nx and A(X) = Yli>j{%i ~ x j)- Note that the spectral deter- 
minant Z N {m) = (— l) Jv nili [(—m 2 ) — x^ is just the characteristic polynomial of 
matrices H from the Laguerre ensemble taken at negative real values of the spectral 
parameter. As is proved in the paper [16] one can express the general correlation 
function of the characteristic polynomials for an arbitrary unitary invariant ensemble 
of P — 2 symmetry class in terms of a (M + L)-sized determinant. The main build- 
ing blocks of that determinant are (monic) orthogonal polynomials vr n (x) = x n + ... 
satisfying 



/ dxw(x)n n (x)ir m (x) = 5 nm cl 

JD 



(5) 



where w(x) is a general weight function, c n are normalization constants and D is 
the corresponding interval of orthogonality. A novel feature revealed in [16] is that 
for M > such a determinant structure contains the Cauchy transforms of the 
orthogonal polynomials 



M e ) 



-/ 

2ni Jd 



w(x) 
ax Tr n (x) 



x — e 



(6) 



alongside with the orthogonal polynomials themselves. For them to be well defined 
we need to have Im{e) ^ 0. 

Actually, the partition function Eq.(3) is given by [16]: 



JC{M f ,M b ) oc 
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A(M b 2 )A(M?) 



(7) 



xdet 
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N-M 



(-k 1} ] 2 ) 
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N-M+l 



(-k 1} ] 2 ) 
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N+L-l 



h N -M (-["^6 M) ] 2 ) h N - M +l (-{lTll M) } 2 ) 

7T N - M (-[m^] 2 ) -kn-m+1 (-[m^] 2 ) 



• • W-it-t^] 2 ) 

• • TTiv+L-i (-[m^] 2 ) 



(-[m 



(i)l2 



TTiV-M+1 (-[mJ. L) ] 2 ) . . . 7T A r +i _i (-[m 



(i)l2 



For the Laguerre ensemble of matrices H with positive eigenvalues X = diag(xi, ...,xn), 
the weight function is just w v (x) = x u e~ Nx , the domain is D = [0 < x < oo] and the 
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monic polynomials are n n (x) = n\ L^(xN), with L v n (xN) being the standard 
Laguerre polynomials. Here v can be taken real valued with v > — 1. To calculate 
the Cauchy transform Eq.(6) we exploit a well-known integral representation for the 
Laguerre polynomials containing the Bessel function J v (x): 

n n (x) = j^j- 2 e Nx x'^ 2 jf dt e-H n+ ^J u (2^Nxl) . (8) 

Let us consider, for definiteness, Im(e) > and further employ the integral repre- 
sentation 



x 



a I" dre- tT(x ~ e) (9) 
Jo 



Then replacing ir n (x) in (6) by (8) and l/(x — e) by (9) we easily perform the inte- 
gration over x first, then integrate over r and arrive at the following representation 
(cf.Eq.(8)): 

h n (e) = ^/^ /2 1 dte~H^H^(2VN^) . (10) 

Here (z) is the Hankel function of the first order. Being actually interested in an- 
alytically continued values of 7i n (x), h n (x) for the region x = —m 2 < we introduce 
the modified Bessel and Macdonald functions according to I v (z) = e~ llTV l 2 J u (iz) 
and K u {z) = fe inu / 2 HW (iz). We then have 

*n(-m 2 ) = J^Jl^ Nm2 ™~ U I dte-H n+ »l 2 I u (2mVM) (11) 
h n {-m 2 ) = [tte-H^ 2 K v (2mVm) (12) 

Substituting such representations into the expression Eq.(7) it is easy to satisfy 
oneself that the right-hand side can be rewritten as (M + L)-fold integral 



JC(Mf,M b ) oc -. -. — p ^3=i l / J 



A(M 6 2 )A(Mf 



detM fe 
det Mr 



(13) 



L L+M 



xf didetitf-MA^e-^RltfWmfNVU)} J[ [t? 2 K v (2mg- L) Ny/ift 

Jti>0 1=1 l=L+l 

Here t > is a diagonal matrix of the size M + L with entries t\, ...,tM+L, and we 
rescaled t — > Nt. Such an equation generalizes the integral representation Eqs.(28)- 
(29) from [15] to nonzero values of v. It is valid for any integer N, L, M. The chiral 
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limit N — > oo can be performed exactly along the same lines as in [15] and the result 
emerging is the one conjectured by Splittorff and Verbaarschot [13], [17]: 

««■ < i4) 

Here Xj = X^ = i : ,^l} and X& = X[i=z,+i,...,L+M} denote the rescaled fermionic and 
bosonic masses respectively as well as Ji — U for i — 1, . . . , L and Ji = Ki for i = 
L+ 1, . . . , L + M. Note that the presence of "compact" (modified Bessel) and "non- 
compact" (Macdonald) functions in the final expressions is a direct consequence 
of the presence of both orthogonal polynomials and their non-polynomial partners 
(Cauchy transforms) in the determinantal representation. One may also wish to 
consider a more general type of potentials V in the probability density, see e.g [7]. 
The related universality questions will be addressed elsewhere [18]. 
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Note that the factor 



is absent in the corresponding expression in [13]. 



" dct X b " 
_dctX f _ 

Instead, it is (implicitly) included in the normalisation of the probability mea- 
sure of the chiral random matrix ensemble. 
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